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Abstract. Specific activator and repressor transcription factors 
which bind to specific regulator DNA sequences, play an important 
role in gene activity control. Interactions between genes coding such 
transcription factors should explain the different stable or sometimes 
oscillatory gene activities characteristic for different tissues. Starting 
with the model P53-MDM2 described into [6] and the process de- 
scribed into [5] we developed a new model of this interaction. Choos- 
ing the delay as a bifurcation parameter we study the direction and 
stability of the bifurcating periodic solutions. Some numerical exam- 
ples are finally given for justifying the theoretical results. 
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1. Introduction 



It is well-known that in cancer, the tumor suppressor protein P53 and its 
partner MDM2 plays an important role. Responses to genotoxic stress are 
regulated by a complex network which center is formed by P53 and MDM2. 
The research effort in P53-MDM2 field is based on the fact that in cancer 
the defects of P53-MDM2 network are almost universal. 

There are 2 types of P53 activity depending of its level: when the level is 
low or there are brief elevation, P53 activates the transcription of the genes 
leading to cell cycle delay at DNA damage-induced checkpoints; when the 
level is high or there are prolonged elevations, p53 activates the transcription 
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of the genes leading to apoptosis. These opposite outcomes are determined, 
at least in part, by cofactors. 

Between p53 and mdm2 exists a reciprocal relationship: p53 stimulates 
the transcription of the mdm2 gene, and thus increase the synthesis of MDM2 
protein; on the other hand MDM2 protein, through its RING domain in the 
C-tcrminus, ubiquinatcs P53 and stimulates its degradation. This mecha- 
nism appears to form a negative feedback loop. We said it appears, because 
recent findings put a question mark on this loop. This arises from the fact 
that there are 2 MDM2 proteins: P90MDM2 and P76MDM2 [7]. P53 acts 
on the promoter P2 of the P76MDM2 and thus stimulates its production but 
P76MDM2 has a truncated N-terminus and can not bind P53 protein and 
thus can not ubiquinate P53, which is needed for degradation of P53. In- 
stead, P90MDM2 has a complete N-terminus so it can bind P53 and thus can 
promote P53 degradation. For the loop to be complete, it requires an interac- 
tion between P76MDM2 and P90MDM2, more precisely P76MDM2 should 
stimulate P90MDM2 binding to P53 or more, P76MDM2 should act on PI 
promoter of P90MDM2 , directly or through an intermediary substance or 
as a cofactor. Because the details of the interactions of the MDM2 isoforms 
(P76 and P90) are only at the beginning to be worked out [7] , and because the 
activators of the PI promoter of the MDM2 (for the P90MDM2 protein) are 
unknown, it could be suspected a link between P76MDM2 and P90MDM2 
that close the loop (in the way presented above). Till now, there are no 
findings to negate this link and moreover, using this loop in our model, we 
obtain an oscillatory behavior similar to that observed experimentally [2,3]. 

We used as a base for our model the model described in [6]. We will 
present that model below: 

yi = i>xi - (Ai + Ai2i/2)?/i 

X2 = Vfi.yiit - t)) - r]2X2 

. m = i'X2 - A22/2 

where: (p is the rate for transcription (maximal) is the rate for translation 
(maximal), r]i,r]2 are the rates for mRNA degradation, Ai, A2, A12, A21 are 
the rates for protein degradation, xi,X2 are mRNA concentrations, yi,y2 
are protein concentrations and Ai is the maximal activation degree for the 
transcription of the p53-gene. 
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Our model is: 



xi = tpAi - r]ixi{t) 

yi = tpxiif?) - (Ai + Ai2?/2(^ - Ti))yi{t) 

m = iJX2{t) - (A2 + X2iyi{t - Tz))y2{t) 

The notations are identical as the previous ones and: ti is the delay required 
for MDM2 to bind P53 plus the time required for the interaction (under 
research) between P76MDM2 - P90MDM2, and also include the time for 
translocation of P53 in cytosol [5] (this is also a mechanism for the down- 
regulation of P53). T2 is the delay required for P53 to enter in the nucleus to 
bind P2 promoter of the mdm2 gene, is the delay required for the HAUSP 
to interact with P53 and MDM2 and to deubiquinate both proteins. A21 is 
degradation rate for MDM2 protein induced by P53. Recent findings show 
that HAUSP (also known as USP7) ,an ubiquitin hydrolase, plays a role in 
P53-MDM2 degradations. Its role, in the presence of P53, is to deubiquinate 
MDM2 and keeps a high MDM2 level. To simplify the expressions that will 
appear in the calculus we use some notations: 771 = 61, Ai = oi, A12 = 012, 
r]2 — b2, X2 — a2, A21 = 021 and also put numerical values for some parameters 
as follows: (p — 1, ip = 1, Ai = 1. These changes have no mathematically 
effect on our system. Finally, we will consider Ti = T2 = T3 = r, the reason 
is that without this hypothesis the calculus become extremely complicated 
and the final result will not differ qualitatively from the calculus with this 
hypothesis. With these specifications made, our system became: 

Xi{t) = 1 — biXi{t), 

yi{t) = xi{t) - (ai + ai2?/2(t - r))yi{t), 

i:2{t) = f{yi{t-T))-b2X2{t), 

hit) = X2{t) - (02 + a2m{t - T))y2{t) 

where / : 1R_|_ — > IR, is the Hill function, given by: 

/(^) = ^ (2) 

with 77. e IN*, a > 0. The parameters of the model are assumed to be positive 
numbers less or equal than 1. 

For Ti = 0, r2 = 0, 021 = in our model, we obtain the model from [6]. 
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The paper is organized as follows. In section 2, local stability for the 
equilibrium state of system (1) is discussed. We investigate the existence 
of the Hopf bifurcation for system (1) using time delay as the bifurcation 
parameter. In section 3, the direction of Hopf bifurcation is analyzed by the 
normal form theory and the center manifold theorem introduced by Hassard 
[4]. Numerical simulations for justifying the theoretical results arc illustrated 
in section 4. Finally, some conclusions are done and further research direc- 
tions are presented. 

2. Local stability and the existence of the Hopf bifurcation. 

For the study of the model (1) we consider the following initial values: 

xi{0) = xi,yi{e) ^ip,{9),9e [-T, 0], X2(0) = X2, y-ziO) = ^2(^), 9 e [-r, 0], 

with (pi,f2 differentiable functions. 
From (1) and (2), we have: 

Proposition 1. If yw is one positive root of the equation: 



+01616200213:^ + (610020162 — 620021)0; — 00262 = 0, 
then the equilibrium state of system (1) is given by: 



(3) 



xio = —, yio, 3^20 = 7-/(2/10), 2/20 = 7-7 : (4) 

61 61 61(02 + 0212/10) 

The proof is immediately. The system should be resolved by nullifying 
the right side of (1). 

We can notice that equation (3) has at least one positive real root. Indeed, 

if: 

k{x) = 0161620210;"+^ + (61O12 - 62021 + ai6i0262)x"'+^ - 0262a;"'+ 
+01616200210;^ + (610020162 — 620021)0; — 00262 = 0, 

then 

k(0) — — 00262 < 0, lim k(x) — 00. 

t — >oo 

So that there is at least one equilibrium point with positive components. 
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We consider the following translation: 

xi^ui + xio, yi^U2 + yio, = % + a^2o, y2^U4 + 2/20- (5) 
With respect to (5), system (1) can be expressed as: 

U2{t) = ui{t) - (fli + ai2y2o)u2{f) ~ ai2yiQUi{t - r) - ai2U2{t)ui{t - r), 
Ui{t) = us{t)-{a2 + a2iyio)u4:{t)-a2iy2oU2{t - T)-a2iU2{t - T)ui{t). 

(6) 

System (6) has (0, 0, 0, 0) as equihbrium point. 

To investigate the local stability of the equilibrium state we linearize 
system (6). We expand it in a Taylor series around the origin and neglect 
the terms of higher order than the first order. Letting X{t) = {yi{t\ V2{t), 
V3{t), V4{t))'^ be the linearized system variables, system (6) can be expressed 
as: 

X{t) = AX{t) + BX{t - t) (7) 

where 
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and pi = /'(2/10). 

The characteristic equation corresponding to system (7) is det{XI — A — 
Be-^'') = which leads to: 

(A + 61) (A^ + h\^ + cX + d+ {gX + h)e-^^^) =0, (8) 

where 

6 = ai + a2 + 62 + 0122/20 + 0212/10, 

C = 62(01 + 02 + O122/2O + O212/I0) + (Ol + Oi22/20)(02 + O212/I0), 

d = 62(01 + Oi22/2o)(o2 + O212/10), 

g = -O122/10O212/20, h = -Oi22/io(&202i2/2o - pi)- 
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First, we consider the case without delay r = 0. According to the Routh- 
Hurwitz criterion we have: 

Proposition 2. When there is no delay, the equilibrium point (xio, yio, 
^20, 2/20) of system (1) is locally asymptotically stable if and only if: 

b>0, {c + g)b>d + h. 

In what follows, we investigate the existence of the Hopf bifurcation for 
system (1), using time delay as the bifurcation parameter. We are looking 
for the values Tc such that the equilibrium point (a;io, yio, 2:20, 2/20) changes 
from local asymptotic stability to instability or vice versa. This is specific 
for the characteristic equation with pure imaginary solutions. Let A = ±iu! 
be these solutions. We assume a; > 0. It is sufficient to look for A = icu, root 
of (8). We obtain: 

{—bu'^ + d + hcos2ujT + ujgsin2ujT — 
—uj^ + cu! + ujgcos2ujT — hsin2ujT — 0. 

By squaring and adding that it follows: 

io"" + {b^ - 2c)lo^ + uj\^ - 2bd - g^) + d^-h^ = Q. (9) 
Let it be ^ = o;^. Then (9) becomes: 

+ hz"^ + hz + /a = 0, (10) 

where 

h = b'^ - 2c, /a = - 2bd - g^, h = d^ - h^. 
Using straight calculus /i > 0, /a > 0. We have: 

Proposition 3. Equation (9) has at least one simple positive root if and 
only if I3 < 0. 

Proof. If Z3 > equation (10) has no real positive roots. If ^3 < then: 

F{z) = / + kz"^ + I2Z + k 

has hm F{z) — 00 and F{0) — I3 < 0. In these conditions equation (9) has 
at least one simple positive root. 

Since delay r is finite, the characteristic equation (8) is a function of 
the delay. Hence the roots of equation (8) are also functions of the delay. 
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Suppose that ujc is the least positive simple root of (8). We will proof that 
there is Tc such that u}{tc) = tOc- At this value of r, the trivial solution will 
lose the stabihty. Defining 

from (8) result 

G{\) = e'^^ (12) 
Thus, we need to find the value of lo such that 

|G(A)| = 1. (13) 

The corresponding values of r are obtained from: 

r = ^[arg{G{iu;)} + 2A;7r], k = 0,1,2,.... (14) 

The critical delay, Tc is the smallest positive value of r satisfying (14). Dif- 
ferentiating (12) implicitly with respect to r, we obtain: 

dX ^ Xg + X^h 

— = 2- 



dr e2^^ (3A2 + 26A - c) - 2r5 + /i - 2t/iA ' 

Then, it is evaluated at A = iuc and r = Tc- 

dX _ ^ ujcgL2 + u^hLi ^ ^. uicgLi + ulhL2 , . 

where 

Li = (— c — 3uj'^)cos2ucTc — 2bUcSin2ucTc — 2gTc + h 
L2 — (— c — 3u!^)sin2uJcTc + 2huJcCos2ujcTc — 2huJcTc- 

Prom the above analysis and the standard Hopf bifurcation theory, we 
have the following result: 

Proposition 4. For r — Tc- 

(dX\ _ uJcgL2+^^hL^ 

and a Hopf bifurcation occurs at the equilibrium state given by (4) as r passes 
through Tc- 
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3. Direction and stability of the Hopf bifurcation 

In this section, we study the direction, stabihty and the period of the 

bifurcating periodic solutions. The used method is based on the normal 
form theory and the center manifold theorem introduced by Hassard [4]. 

From the previous section, we know that if r = Tc then any root of (8) of 
the form 

A(t) = Q!(t) ± iuj{T) 

satisfies a{Tc) = 0,u{tc) = ujq and \t=tc 0. 

dr 

Define the space of continuous real- valued functions as C = C([— Tc, 0], H^). 
Expanding / in (6) in Taylor series around (0, 0, 0, 0)^ we can rewrite system 
(6) in the form: 

X{t) = AX{t) + BX{t - r) + F{X{t),X{t - r)) (16) 

where 

F(X(t),X(t - T))=iO,F\u2it)Mt - r)\F\u2{t - T)\F\u2{t - T\u^{t))f , 

F'^{u2it),Ui{t - T))=-ai2U2it)u4{t - t), 

F^{u2it - T))=|p2«2(t - + |p3«i(t - r), 

F\u2{t - T),U4{t)) = -a2lU2{t - T)u4t),p2 = f"{yio),P3 = fiVlo)- 

In r = Tc + e M, we regard fj, as the bifurcation parameter. For 
$ e C we define a linear operator: 

L(//)$ = ^$(0) + 5$(-r) 

and a nonlinear operator: 

= (0,F2($2(0),$4(-r)),F3($2(-r)),F^($2(-r),$4(0))r. (17) 

For $ e Ci([-Tc,0],(D*^) we define: 

f d<l>{9) 



de ' 

[ A$(0) + 5$(-Tc), ^ = 0, 
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0, ee[-Tc,0) 

and for ^ e ^"^([0, rJjCD*^), we define the adjoint operator A* of A by: 



A*^{s) 



Then, we can rewrite (16) in the following vector form: 

Xt = A{fx)Xt + R{i^)Xt 

where Xt = X{t + 9) for 6 e [-Tc,0]. For $ e ( [-Tc, 0],(D*^) and * e 
C([0, rc],(D*^) we define the following bilinear form: 



< 



^(^),$(^) >= ^(0)$(0) + r ^(^ + rc)B$(^)(i^. 



Then, it can verified that A* and ^(0) are adjoint operators with respect 
to this bilinear form. 

In the light of the obtained results in the last section, we assume that 
±iuc are eigenvalues of A{0). Thus, they are also eigenvalues of A*. We can 
easily obtain: 

^e) = ve^'^, ee[-Tc,o] (18) 

where v = (vi, V2, vs, v^f , 

^ Ai + 62 Air 1 (Al + 62)(Ai +01 + 0122/20) 2Air 

Pi Piauyio 
is the eigenvector of ^(0) corresponding to Ai = iu>c and 

*(s) = we^^', se[0,rc] 

where w = {wi, W2,W3, W4), 

1 + Ai ai2yioe-^i^(fei + Ai) 

V r] (02 + 021^/10 + Ai) (62 + Ai)77 

ai2ywe-^^^{bi + Ai) 

(02 + 0212/10 + Ai)?7 



/, , , V a^yioe ^^^{h + Ai) auVwe ^^^{bi + Ai) 
V = -^2(01 + Ai - ■ . . w, . . . vs — — V4+ 

(a2 + a,2iyw + Ai)(62 + Ai) 02 + Ai + a2i?/io 

(02 + 021^10 + Ai) (02 + Ai) 02 + Ai + a2iyio 

is the eigenvector of A* corresponding to A2 = —iuJc- 

We can verify that: < ^(s),$(s) >= 1, < ^(s),$(s) >=< ^(s),$(s) >= 
0, < >= 1. 

Using the approach of Hassard [4], we next compute the coordinates to 
describe the center manifold ji — 0. Let Xf — X{t + 9), 9 e [— Tc, 0], be 

the solution of equation (7) when /i — and 

z(t) =< ^, Xt >, w{t, 9) = Xt{e) - 2Re{z{t)<!?{9)}. 

On the center manifold fig, we have: 

w{t,9) = w{z{t),z{t),9) 

where 

W{Z, Z, 9) = W2o{9)— + WiiZZ + Wo2{9)— + W3o(6')— + . . . 

in which z and z are local coordinates for the center manifold in the 
direction of ^ and ^ . 

For solution Xt of equation (7), since = 0, we have: 

z{t) ^ \iz{t) + g{z,z) (19) 

where 



g{z, z) = *(0)F(w(z, z, 0) + i?e(z(t)$(0))) 

2p' Z^Z 

= 920— + giizz + 5o2y + 921— + 



(20) 



From (17), (18) and (20) we obtain: 

920 = F^qW2 + F2^o^3 + F^QWi, gu = Ff^ ?/;2 + Ff^w^ + F^^W4, 

902 = F^^W2 + F^^W3 + F(|W4, 



10 



where 



f2 

-^02 

-^11 
-'^20 
F02 



P2vle-^^'^, 



P2V2e-^^2r 



P2^^2^^2, F62 

-2a2i'i;2'J^4e~'*'^'^, = -a2iV2V^e~^^'^ + ■i;4t;2e~'^2'^, 
-2ai2^;2'J^4e~'*'2'^, 



and 



where 



921 = i^2i^2 + -^21^3 + ^21^4 



(22) 



-ai2(2t;2wfi(-r) + V2W^q{-t) + t;4«;io(0)e"^2^ + 2vAwli{-T)e-^^^) 
F|i = p2(2^J2«'?i(-^)e"^'^ + V2wlQ{-T)e-^^^) + p3t^2^2e-^i^ 

= -a2i(2v2«;fi(-r)e-^i" + ^2tx;^o(0)e-^^" + t^4W^2o(-r) + 2v4U)?i(-r)) 



21 
with 



<i-r) = -^^2e-^- - |^^2e-^- + ^2^6"^^- 

^2^(_r) = ^^26-^^^ - ^vle-'^^ + El 

-^^;4e-^- - |f^^)4e-^- + i?2'e-2^- 

t(;4^(-r) = ^^46-^1^ - ^^v^e-^^^ + 
Ai _ Ai 

2 /n\ ^20 902 - . rp2 

«^20(0) = -^^2 - ^^^2 + ^2 

2 /^^ ^11- I 772 

^11 (0) = -^-^2 - -— ^^2 + E^ 

M Al 

4 /^x 5*20 5'20 - , Tp4: 

«^2o(0) = - ^f4 + ^2 

4 /n^ •S'll ^11- I 774 

<i(0) = -;-^4 - ^^^4 + E^, 
Ai Ai 

Ef^ Ef respectively are the components of the vectors: 



Eo 



-{A + e-2^i^=S - 2\iI)-^F: 
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E,^-{A + B)-'Fu, 



where F20 = (0, F|o, F|o, Fi,f, Fn = (0, Ff,, Ff„ Ff, 



'4 \T 
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Based on the above analysis and calculation, we can see that each gij in 
(21), (22) is determined by the parameters and delay from system (1). Thus, 
we can explicitly compute the following quantities: 



= - ReX.in) '^^ = '^^ = 2i?e(C,(0)). 

(23) 

In summary, this leads to the following result: 

Proposition 5. In formulas (23), 112 determines the direction of the 
Hopf bifurcation: if > 0(< 0), then the Hopf bifurcation is supercritical 
(subcritical) and the bifurcating periodic solutions exit for r > Tc{< Tc); /?2 

determines the stability of the bifurcating periodic solutions: the solutions are 
orbitally stable (unstable) if P2 < 0(> 0); and T2 determines the period of the 
bifurcating periodic solutions: the periodic increase (decrease) ifT2 > 0(< 0). 

4. Numerical examples. 

For the numerical simulations we use Maple 9.5. In this section, we 
consider system (1) with ai = 02 = 0.13, = 0.02, 021 = 0.02, bi = 0.8, 
62 = 0.01, a = 4; ai2 = 021 because there is molecular interaction between 
MDM2 and P53, one molecule to one molecule. Waveform plot and phase 
plot are obtained by the formula: 

X(t+^) = z(t)*(^)+z(t)i>(^)Hi«;2o(^)-z'(t)+wii(^)z(t)^(t)H^wo2(^)^(i)'+Xo, 

where z{t) is the solution of (19), $(6*) is given by (18), W2o{0) , wii{6) , Wo2{0) 
are given by (22) and Xq — (xio, ?/io, 3:20; ?/2o)'^ is the equilibrium state. 

Forn = 2 we obtain:xio = 1.2500000, |/io = 0.72279716, 2/20 = 79.96962531, 
X2Q = 11.55208766, /i2 = -15.56012572, /?2 = -0.00020024, T2 = -0.169703418, 
a; = 0.01173958, r = 90.21567180. Then the Hopf bifurcation is subcritical, 
the solutions are orbitally stable and the period of the solution is increasing. 
The wave plots and the phase plot are: 
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Waveplot (t, yi) 
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Phaseplot (^1,2/2) 
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Forn = 4 we obtainixio = 1.2500000, ?/io = 0.82091152, ?/2o = 69.63487984, 
X2o = 10.19581588, //2 = -22.21740930, /?2 = -0.00987558, r2 = -0.86252133, 
u = 0.02969208, r = 26.61818721. Then the Hopf bifurcation is subcritical, 
the solutions are orbitally stable and the period of the solution is increasing. 
The wave plots and the phase plot are: 



Waveplot {t,yi) 



Waveplot {t,y2) 



Phaseplot (^1,^2) 






200 400 600 



For n = 163 we obtain:xio=l. 2500000, yio=0.99390609, 2/20=56.38320475, 
X20 = 8.45060883, fi^ = -12.63855144, P2 = -0.53197047, T2 = 7.14847952, 
llj = 0.42317766, r = 0.00213625. Then the Hopf bifurcation is subcritical, 
the solutions are orbitally stable and the period of the solution is increasing. 
The wave plots and the phase plot are: 
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For n = 164 we obtam:xio=l. 2500000, 2/io=0.99394289, 2/20=56.38087608, 
X20 = 8.45030131, /i2 = 4.70953378, /?2 = -1.32779287, = 3.42681695, 
cu = 0.42448028, r = 7.40096599. Then the Hopf bifurcation is supercritical, 
the solutions are orbitally stable and the period of the solution is increasing. 
The wave plots and the phase plot are: 
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Phaseplot (2/1,^2) 
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4. Conclusions. 

Recent dynamic studies of P53 and MDM2 proteins suggest that their 
responses in individual cells have cyclic behavior and their characteristics 
are compatible with a digital clock [2]. Similar behavior we obtained in our 
mathematical model. The qualitative result confirms the existence of a limit 
cycle which can be characterized using the coefficients from (23). 

The model described by (1) suggest the investigation of the model defined 
by the following system of nonlinear differential equation with distributed 
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delay: 

Xi{t) = 1 — biXi{t), 

yiit) = xi{t) - (ai + ai2?/2(t))?/i(t), 
1 

i2{t) = - / f{yi{t - s))ds - b2X2{t), 
y2{t) = X2{t) - (02 + a2iyi{t))y2{t) . 
The model will be studied in a future paper using the method from [1] . 
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